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ABSTRACT 


This  report  examines  gas  dynamic  stabilization  of  high  pressure, 
highly  non-uniform  electrical  discharges  for  two  electrode  geometries. 

The  geometries  studied  are  the  multiple-point  to  plane  and  the  multiple- 
wire  to  plane.  A  successful  analytical  model  is  devt  loped  utilizing  the 
anode  extension  concept,  which  describes  the  interactions  and  stabilizing 
effects  of  flow,  separating  the  effects  of  convection  and  turbulence,  in 
the  point-plane  geometry.  Comparisons  of  theory  and  data  are  presented, 
illustrating  actual  observed  and  calculated  breakdown  voltages  and 
currents  for  a  number  of  cases.  Flow  stabilized  high  power  density 
electric  discharges  are  applicable  to  high  power  electric  discharge  lasers, 
magnetohydrodynamic  power  generation,  and  other  new  systems  and  technolo¬ 
gies. 


vi  ii 


I.  INTRODUCTION 


A.  Background 

Electric  discharges  have  been  the  subject  of  intensive  study  for  many 
years.  Numerous  books,  papers,  reports,  etc.  have  been  published  which 
characterize  and  analyse  the  various  complex  and  interrelateu  processes 
which  take  place  in  discharges  of  several  gas  compositions,  pressures,  and 
electrode  geometries.  The  primary  objective  in  most  of  these  studies  has 
been  to  attempt  to  relate  the  atomic  level  processes  to  the  observed 
macroscopic  discharge  properties  such  as  breakdown  voltage,  current  and 
overall  discharge  stability  for  use  in  engineering  analysis  and  application 
in  design.  Some  of  the  more  established  applications  of  this  type  of 
gaseous  electronics  include  glow  discharge  lighting,  the  high  voltage  arc, 
corona  discharge  analysis  around  power  lines  and  power  line  design,  high 
voltage  switching,  spark  plug  and  spark  gap  design,  thermal  printing 
devices  and  many  others. 

More  recent  applications  of  this  electric  discharge  technology  have 
been  generally  aimed  at  systems  which  are  geared  to  creating  very  high  power 
densities  within  the  discharge,  thrusting  physicists  and  engineers  headlong 
into  a  direct  confrontation  with  the  stability  problem,  not  just  with  the 
idea  of  "defining  the  envelope",  to  use  an  aeronautical  term,  but  to  develop 
ways  of  expanding  the  natural  discharge  stability  boundaries  through  any 
practical  means.  The  recent  high  power  applications  of  which  we  speak 
include  magnetohydrodynamic  (MHD)  power  generation,  plasma  chemistry, 
electric  discharge  convection  lasers.  E-beam  ionizer-sustainer  lasers  and 


others. 


As  a  result  of  this  growing  requirement  for  increasingly  higher  power 
densities,  many  analyses  have  been  done  in  the  stability  area  and  much 
insight  has  been  gained  in  the  many  modes  and  regimes  from  which  discharge 
instability  can  develop  to  cause  a  severe  local  breakdown  in  a  delicate 
balance  of  parameters.  One  of  the  successful  techniques  which  has  appeared 
in  the  engineering  of  high  power  density  systems  is  the  pulsed  system 
where-in  the  pulse  time  is  shorter  than  the  characteristic  times  of  the 
various  instability  modes  which  must  be  avoided.  The  E-beam  and  pulsed 
ionizer-sustainer  lasers  are  examples  of  the  success  of  this  approach 
(Reference  1).  However,  present  E-beam  systems  are  large,  expensive  and 
delicate.  In  addition,  in  order  to  achieve  high  average  power  output 
from  the  pulsed  systems,  while  the  systems  can  be  scaled  up  in  size  to  any 
practical  limit,  once  the  size  limit  has  been  reached,  the  only  other 
alternative  is  to  increase  the  pulse  repetition  frequency  (PRF) .  But  this 
option  is  also  limited,  usually  by  a  characteristic  optical  cavity  clearing 
time  which  depends  on  cycling  out  by  flow  the  products  of  one  pulse  before 
initiating  the  next  one.  Failure  to  allow  sufficient  clearing  time  will 
result  in  cumulative  space  charge  and  gas  heating  instabilities  with 
resultant  low  power  breakdown. 

Due  to  the  inherent  limitations  cited  above  and  the  implicit  depen¬ 
dence  of  those  systems  on  gasdynamics,  the  present  research  has  been  pur¬ 
sued  in  order  to  better  define  the  role  of  turbulence  in  a  high  pressure 
discharge  as  has  been  observed  by  many  researchers  (References  Z,  3,  *♦,  5 
and  6).  An  understanding  of  the  fundamental  modes  of  gasdynami c  coupling 
which  have  been  observed  to  affect  breakdown  voltage  and  current  may  lead 


to  new  ways  to  improve  present  high  power  density  systems  through  bulk  or 
electrode  boundary  layer  stabilization  techniques;  and  may  also  point  the 
way  for  development  of  sufficiently  high  power  self-sustained  CW  electric 
discharge  systems  which  are  inherently  more  compact,  rugged,  and  less  ex¬ 
pensive  than  the  types  mentioned  above. 

B.  Present  Work 

Section  A  of  Part  II  (Fundamentals  of  Gasdynamic  Effects  on  a  Highly 
Non-Uniform  Electrical  Discharge)  of  this  report  develops  in  detail  a 
theoretical  model  for  the  calculation  of  breakdown  voltage  of  a  *i- 
point  to  plane  (points  positive)  electric  discharge  based  on  the  '‘anode 
extension"  concept  initially  proposed  in  Reference  7  and  further  discussed 
in  Reference  3.  The  model  incorporates  flow  and  turbulence  effects  and  is 
similar  in  physical  basis  to  the  "corona  cloud"  model  of  Hepworth,  Klene, 
and  Tozer  discussed  in  Reference  8  for  a  static  (non-flowing)  non-uniform 
electric  discharge.  The  present  model  successfully  separates  the  gasdy- 
namics  into  the  pure  convective  and  pure  turbulence  effects,  and  isolates 
specific  electrode  geometry  and  gas  parameters  in  individual  terms.  This 
facilitates  understanding  of  the  fundamental  relationships  between  the  gas 
properties,  electrode  orientation  and  geometry,  and  the  individual  gasdy¬ 
namic  coupling  mechanisms. 

This  section  also  presents  a  qualitative  discussion  of  the  breakdown 
characteristics  of  the  uniform  and  non-uniform  discharge  geometries  high¬ 
lighting  the  similarities  in  a  qualitative  and  semi -quant i tat i ve  manner. 

It  is  shown  that  Meek's  breakdown  criterion  for  the  non-flowing,  uniform 
field  and  the  equation  developed  in  Section  A  for  the  non-uniform,  flowing 
discharge  transform  into  very  nearly  identical  expressions  in  the  middle 
ground  of  a  flowing  uniform  field.  This  illustrates  the  general  nature  of 
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A 

Che  equations  of  Section  A  and  provides  straightforward  physical  insight 
into  the  final  breakdown  mechanisms  in  the  geometric  extremes.  Radial 
diffusion  and  spreading  is  considered  as  related  to  the  mean  or  "drift" 
field  in  a  highly  non-u\iform  geometry. 

Section  B  of  the  report  discusses  the  typical  discharge  section  used 
in  Section  A  with  the  three  rows  of  pins  replaced  by  three  thin  wires  (see 
Fig.  lb).  A  preliminary  analysis  which  was  done  in  the  terms  of  Reference 
9  basically  indicated  the  potential  for  improving  discharge  input  power  by 
a  factor  of  2-3  with  the  wire  anode.  Various  degrees  of  discharge  input 
power  enhancement  have  been  noted  in  similar  configuration  modifications 
in  connection  with  MHD  laser/power  generation  applications  (Reference  10). 

Recent  experimental  results  in  this  area  are  presented  and  discussed 
in  relation  to  the  points-plane  configuration.  To  date  due  to  power 
supply  limitations  no  improvement  in  performance  has  been  noted,  and  some 
recommendations  for  further  research  in  this  area  are  presented,  when  a 
larger  power  supply  is  available. 

Section  C  of  Part  II  discusses  electrode  design  geometry  and  electric 
field  distributions  for  optimum  electron  production  and  for  best  reduced 
electric  field  strength  (E/N  or  E/P)  for  laser  pumping. 

Part  III  presents  a  data  summary  of  current  and  past  point-plane  and 
wire-plane  laboratory  work,  reviews  the  theory  and  experimental  trends, 
indicating  important  correlations  and  features  of  the  anode  extension  model. 
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II.  FUNDAMENTALS  OF  GAS  DYNAMIC  EFFECTS  ON  A  HIGHLY 
NON-UNIFORM  ELECTRICAL  DISCHARGE 


A.  Prediction  of  Flow  and  Turbulence  Effects  on  Breakdown  Voltage 
1 .  Description  of  Discharge  Geometry  for  Analytical  Formulation 

A  significant  body  of  breakdown  data  has  been  generated  in  var¬ 
ious  highly  non-uniform  electric  field  geometries  with  positive  and  neg¬ 
ative  polarities  at  the  highly  stressed  electrode.  The  point-to-plane 
geometry  with  the  points  positive,  however,  has  been  the  particular  case 
through  which  the  present  analytical  model  has  been  developed;  but  the 
model  is  sufficiently  general  that  the  basic  flow  and  electric  field 
interactions  described  should  also  be  useful  in  analysis,  interpreta¬ 
tion,  and  trouble  shooting  of  many  other  non-uniform  geometries  and 
systems  both  with  and  without  flow. 

The  specific  apparatus  used  in  the  current  research,  described 
in  detail  in  Reference  11  consists  of  multi-point  to  plane  electrode 
arrangement  (see  Fig.  la)  which  is  mounted  in  a  wind  tunnel  with  avail¬ 
able  test  section  velocities  on  the  order  of  200  m/s  of  ambient  pressure 
air.  Varying  degrees  of  turbulence  can  be  introduced  by  placing  turbu¬ 
lence  generating  screens  upstream  of  the  electrode  set-up.  Two  flow 
modes  are  discussed  ir  this  report,  namely, 

a.  Normal  Flow  -  directed  through  the  turbulence  screen,  past 
the  positive  multi-pin  anode,  through  the  interelectrode  gap,  and  thence 
through  the  grounded  cathode,  which  is  designed  to  approximate  a  plane 
with  low  blockage  to  permit  maximum  gas  flow  rates  (see  Fig.  2a). 

b.  Reverse  Flow  -  same  as  above  except  that  the  anode  and  cat¬ 
hode  positions  are  reversed  (see  Fig.  2b)  and  the  turbulence  screen  was 
omi tted. 
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While  our  earlier  effort  which  included  cross-flow  tests  (Fig. 
2c)  was  based  on  an  assortment  of  turbulence  screens  which  produced  a 
variety  of  spectra  and  intensities,  the  present  work  (which  includes 
verification  of  past  trends,  and  additional  data  taking  for  better  curve 
definition  at  low  speeds)  is  done  utilizing  only  screen  #9  design.  This 
screen  proved  to  be  most  successful  In  the  prior  work  (References  3,  7, 


2.  General  Description  of  the  Discharge 

With  no  voltage  applied  across  the  gap,  and  no  gas  flow,  there 
exist  a  small  number  of  random  electrons  present  in  the  gap  which  are 
produced  by  outside  sources  (Reference  12).  As  voltage  is  applied  these 
electrons  under  the  influence  of  the  resultant  electric  field  begin  to 
move  toward  the  anode  producing  a  small  saturation  current  whose  magni¬ 
tude  is  limited  by  the  number  of  electrons  present  and  their  average 
drift  velocity  in  the  field.  In  a  highly  asymmetric  geometry  such  as 
the  point-to-plane,  the  electric  field  produced  is  very  non-uniform, 
being  highest  at  the  point  surface  and  falling  off  rapidly  in  a  hyper¬ 
bolic  fashion  with  distance  from  the  point  as  illustrated  by  the  axial 
field  distribution  in  a  hyperboloi dal  point-to-plane  gap  (Reference  17) 
which  is  assumed  in  the  following  analysis: 


EW  -^Ep 

R  =  point  radius 
x  =  axial  distance  from  point 


Ep  =  field  strength  at  point  = 
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V  =  applied  voltage 
d  =  gap  length 

Oue  to  this  field  asymmetry,  as  the  voltage  increases  further,  the  field 
strength  at  the  point  soon  reaches  a  value  which  permits  the  free  elec* 
trons  moving  toward  the  point  to  gain  enough  energy  to  ionize  by  colli¬ 
sion  some  of  the  gas  atoms  and  molecules  near  the  point  surface,  thus 
additional  electrons  and  positive  ions  are  produced  resulting  in  a  rise 
in  the  current.  The  voltage  at  which  these  first  few  ionization  events 
occur  is  called  the  corona  initiation  voltage  and  is  usually  designated 

V  ,  i.e.,  when  Ep  Z  E  ,  Vs s  V  : 

c  c  c 


2V 

Ep  ■*  Ec  ~  Rlog(4d/R) 


(2) 


EcRlog(Ud/R) 


(3) 


From  Reference  13  the  threshold  field  Ec  can  be  calculated  from  Peek's 
Law  (1929)  modified  by  Gary  et  al  (1972)  for  a  poi nt-to-pl ane  type 
geometry,  i.e., 


Ec  -  Eom 


I  + 


K 


(a|) 


where 


w 


Eq  -  31  KV/cm  (near  spark  breakdown  of  air) 

M  =  coefficient  describing  conductor  surface  (0. 6-1.0) 
K  =  0.308  cm* 

R  =  radius  of  point 

A  =  0.392  P/(273+T)  =  1  for  air  at  STP 
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P  2  pressure  in  torr 

_  -0r 
T  =  temperature  tn  C 

Continued  Increases  in  voltage  moves  the  high  fie  id  region  further  and 
further  into  the  gap.  At  the  position  where  the  reduced  field  value 
(E/P)  is  approximately  equal  to  that  required  for  ionization  by  elec¬ 
tron  collision  (Ej/P)  ,  free  electrons  will  begin  to  produce  avalanches 
into  the  positive  point.  From  this  critical  distance  x£  ,  each 
avalanche  will  produce  an  exponential  increase  in  electrons  and  positive 
ions  left  in  the  aftermath.  The  number  of  electrons  produced  in  the  non- 
uniform  field  is  given  by  the  well  known  expression  (Reference  12) 

x2 

/x  adx 

n  -  nQe  1  (5) 

where 

nQ  2  original  number  of  electrons 

a  5  Townsend’s  first  ionization  coefficient 
x  2  avalanche  travel  distance 

The  coefficient  a  is  a  strong  function  of  the  electric  field  as  illus¬ 
trated  by  the  usual  semi -analyti cal  expression  found  in  many  references  - 

BP 

a  .  APeE^x)  (6) 

where  A  and  B  are  gas  constants,  p  is  pressure.  It  is  necessary, 
therefore,  to  integrate  the  exponent  since  in  the  region  under 


discussion  E  and  therefore  a  vary  greatly  with  x  as  the  avalanche 
rushes  toward  the  point,  Fig.  3*  An  average  or  effective  ionization 
coefficient  can  be  defined  as  follows: 

x 

“d*  (7i 
o 

X 

o 

which  will  produce  the  same  number  of  particles  in  traveling  the  same 
distance,  and  then  we  can  write 

a(x-xQ) 

n  *  n  e 
o 

similar  to  the  simpler  uniform  field  case.  In  the  following  discussion 
the  bar  will  be  dropped  for  convenience,  and  xq  assuned  to  be  0.  Due  to 
the  much  greater  electron  drift  velocities,  the  electrons  from  these 
avalanches  are  quickly  drawn  into  the  positive  point  leaving  the  positive 
ions  behind  in  a  maelstrom  of  activity  including  electrostatic  repulsion 
from  the  positive  point,  diffusion  in  all  directions,  neutralization  by 
recombination  with  other  electrons  or  charge  transfer  from  negative  ions, 
etc.  All  of  these  processes  tend  to  dissipate  the  space  charge  concen¬ 
trations  and  smooth  the  non-uniformity  of  the  electric  field  which  is 
generated  by  the  dense  space  charge  near  the  point,  and  by  the  point 
i tsel f. 

However,  if  the  applied  voltage  is  raised  high  enough  (near 
breakdown),  these  dissipative  processes  lag  behind  production  processes, 


and  the  space  charge  density  and  field  intensity  can  become  of  the  same 
order  of  magnitude  as  the  applied  field.  When  this  happens,  the  portion 
of  the  gap  external  to  the  high  field  region  is  no  longer  able  to  dis¬ 
tinguish  the  applied  field  nor  the  position  of  the  point,  and  the  dense 
space  charge  actually  assumes  the  role  of  the  point  forming  a  "virtual" 
anode  (Reference  14)  effectively  diminishing  the  gap  as  seen  from  the 
cathode.  The  effect  of  this  "extension"  of  the  point  is  to  require 
dissipation  of  the  applied  voltage  over  a  shorter  gap,  thus  raising  the 
"drift"  (low  field)  portion  of  the  gap  to  a  uniformly  higher  level  (Fig. 

4).  Reference  14  illustrates  how,  in  a  non-uniform  geometry,  the  presence 
of  space  charge  can  modify  an  originally  hyperbolic  field  distribution 
to  produce  a  uniform  field  over  most  of  the  gap  in  concentric  cylinders 
where  the  inner  electrode  radius  is  much  less  than  the  outer. 

Another  important  effect  which  must  be  accounted  for  is  the 
radial  diffusion  of  the  space  charge  in  the  ionization  region  which 
essentially  distorts  the  shape  of  the  original  point,  making  the  resultant 
virtual  anode  appear  extended  in  the  radial  as  well  as  the  axial  dimen¬ 
sion  (see  Fig.  5) . 

It  will  be  seen  in  the  following  analysis,  that  accounting  for 
the  two  effects  (shortened  gap,  and  modified  radius)  provides  much  in¬ 
sight  into  various  discharge  phenomena  heretofore  not  well  explained 
and  leads  to  a  very  general  form  of  the  breakdown  equation  for  a  highly 
non-uniform  geometry  In  a  turbulent  flow  field.  The  breakdown  equation 
Is  derived  by  accounting  for  modifications  by  the  flow  of  the  quasi- 
statlc  space  charge  distributions  which  develop  in  the  non-flowing  run¬ 
ning  discharge  case.  The  trends  revealed  and  results  obtained  are 
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fundamental  and  adaptable  to  different  gas  mixtures,  geometries,  and  flow 
field  parameters. 

3.  Analytical  Development  of  the  Breakdown  Equation  with  Flow 
Interactions 

The  exponential  charged  particle  multiplication  process  in  high 
reduced  fields  has  already  been  mentioned.  Since  the  current  is  produced 
by  the  motion  of  these  particles  between  the  electrodes,  the  current 
density  can  be  written  as  a  summation  of  the  individual  species  contribu¬ 
tions.  For  the  present  it  is  assumed  that  the  positive  ion  and  electron 
currents  dominate  the  discharge  (Reference  8).  The  current  density  will 
then  be  given  by 

J  "  J|  "  Je  *  enivi  +  eneve  (8) 


where 

j  S  local  current  density 

.  ]  q 

e  =  electronic  charge  (1.6x10  coulombs) 
v.  =  ion/electron  drift  velocities 

I  *6 

n.  =  local  charge  densities 

I  ,6 

Considering  only  the  high  field  ionization  region,  if  each  elec¬ 
tron  entering  causes  an  electron  avalanche,  the  average  number  of  elec¬ 
trons  produced  Is  equal  to  the  multiplication  factor  plus  one,  while  the 
number  of  positive  ions  produced  would  be  just  the  multiplication  factor. 
The  multiplication  factor  is  usually  so  high,  that  the  number  of  positive 
ions  and  electrons  created  in  this  region  are  Identical.  Therefore 
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n.  =  ne  =  n  (9) 

Due  to  the  relatively  low  mobility  of  the  positive  ions,  the 
majority  of  electrons  continue  into  the  anode  and  a  region  of  growing 
positive  space  charge  is  produced  in  the  immediate  vicinity  which  will 
slowly  begin  to  drift  away  due  to  electrostatic  repulsion  from  the  anode. 
Under  these  conditions  the  current  density  can  be  estimated  by  the 
following 

j  *  en  (vj+vg)  (10) 

Under  conditions  of  low  positive  space  charge,  vg  is  much  greater  than 
v.  and  the  equation  reduces  to 

j  =  enve 

But  if  the  positive  space  charge  field  becomes  of  the  same  order  of  mag¬ 
nitude  as  the  local  applied  field,  electrons  may  be  retarded  or  even 
trapped  by  the  space  charge  cloud.  In  this  case  the  electron  drift 
velocities  may  be  reduced  to  the  same  order  of  magnitude  as  the  ion 
velocities,  or  to  zero  and  even  negative  if  the  electron  is  trapped  and 
dragged  in  the  opposite  direction  by  high  positive  ion  fields.  In  this 
case,  the  current  density  approaches 

j  »  env.  (11) 

and  the  particle  density  may  then  vary  between 
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T 


n 


to  n 


do 


It  is  now  proper  to  return  to  the  "anode  extension"  concept  introduced 
above.  As  indicated,  when  the  space  charge  field  near  the  anode 
approaches  the  same  order  of  magnitude  as  the  applied  field,  flux  lines 
from  charges  within  the  remainder  of  the  gap  and  from  the  cathode  will 
no  longer  be  able  to  reach  the  actual  point,  but  will  instead  terminate 
on  charges  in  front  of  the  point.  This  essentially  shortens  the  gap  and 
to  a  first  approximation,  results  in  dissipation  of  the  applied  potential 
over  a  reduced  distance,  raising  the  mean  electric  field  over  the 
remainder  of  the  gap  (Fig.  5). 

In  order  to  continue  the  development,  an  estimate  of  the  length 
of  the  anode  extension  region  (hereafter  designated  by  a  6  )  is 
required.  A  reasonable  estimate  of  the  beginning  of  the  region  of  dense 
positive  charge  would  be  the  point  at  which  the  applied  electric  field 
approaches  the  breakdown  value  for  the  type  of  gas  mixture  involved.  In 
our  example  of  atmospheric  air,  it  is  assumed  that  the  anode  extension 
in  a  non-flowing  discharge  begins  at  the  value  E(x)  *  25,000  volts/cm 
(E/P  s  35).  This  assures  that  the  majority  of  electrons  approaching 
this  position  avalanche  toward  the  positive  point,  depositing  the  slower 
moving  positive  ions  to  form  the  anode  extension  (see  Fig.  5). 

As  soon  as  a  space  charge  concentration  develops,  however,  there 
are  numerous  forces  and  processes  which  act  on  and  within  it  to  dissipate 
it.  Among  the  important  process  to  be  considered  are: 

(a)  diffusion  in  all  directions, 
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(b)  electrostatic  repulsion  by  the  electrode  of  like  sign,  in 


this  case  the  nearby  point, 

(c)  attraction  by  the  opposite  electrode  which  may  be 
insignificant  in  this  region  in  very  non-uniform  fields, 

(d)  recombination  with  supplementary  electrons  supplied  from 
various  sources,  i.e.,  photoelectrons,  cathode  emission, 
free  electrons,  detachment  from  negative  ions,  etc. 

In  a  "running"  discharge,  i.e.,  one  which  has  reached  a  pseudo-steady 
mode  of  operation,  a  quasi-equilibrium  condition  is  ultimately  reached 
in  which  the  charge  production  and  dissipation  mechanisms  locally  balance 
each  other.  The  result  of  this  balance  is  a  quasi -steady  anode  extension 
with  a  finite  and  calculable  extent;  a  distorted  radial  dimension  of  the 
virtual  anode  due  to  radial  diffusion,  and  an  increased  mean  electric 
field  in  the  drift  region  of  the  discharge  (see  Fig.  5).  Reference  3 
alluded  to  these  two  important  dimensional  changes. 

To  calculate  the  equilibrium  boundary  of  5  in  the  axial 
dimension  between  anode  and  cathode,  electrostatic  repulsion,  convection 
(in  the  flowing  case)  and  recombination  must  be  balanced,  diffusion  being 
small  in  the  high  field  region  compared  to  the  axial  component  of  ion 
drift.  The  radial  boundary  is  determined  by  radial  diffusion,  and  re¬ 
combination,  and  in  the  case  of  crowded  multi-pin  electrodes,  by  symmetry 
considerations. 

a.  The  Anode  Extension  -  Axial  Dimension  -  6 

Since  the  length  of  6  is  determined  primarily  by  a  balance 
between  the  ion  drift  out  of  the  dense  space  charge  region,  and  the  dis¬ 
sipation  of  that  charge  by  recombination,  an  expression  for  recombination 


is  required.  The  form  assu»cJ  here  as  in  Reference  3  is  the  three-body 
solution  of  Hinnov  and  Hirschberg,  namely  that  the  recombination  time 


t 


REC 


02) 


whe  re 


P 

n 


'■REC 


=  recombination  coefficient 
=  particle  density 
=  characteristic  recombination 


time 


Now,  the  positive  space  charge  created  in  the  vicinity  of  the 
anode  will  be  subject  to  electrostatic  repulsion  until  it  recombines  with 
an  electron.  Therefore,  the  travel  distance  of  an  individual  ion  during 
its  1 i fetime  will  be 


x 


d 


vi  tREC 


(13) 


where  v.  H  k.E  =  ion  drift  velocity.  Inside  this  drift  distance  will 
be  much  positive  space  charge,  and  outside  it,  most  of  the  charge  will 
have  been  neutralized  by  recombination.  Therefore,  since 


xd  =  6 


(14) 


combining  Equations  12,  13,  and  >4 


S  =  v,  (-Ly) 

2pn 

Substituting  for  n  from  Equation  11 

v.  ev„  2  v.  ev,  2 

5  =  *<-r>  t0 


(15) 


(16) 
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So  far,  flow  has  not  been  introduced  so  the  above  equation  is 
general.  Taking  flow  into  account,  the  following  notation  will  be  used: 

=  anode  extension  in  flowing  gas 

5  =  anode  extension  with  no  flow 

o 

6 

6*  =  flow  value  normalized  by  no  flow  value  (j— •) 

o 

Typical  flow  velocities  in  practice  range  up  to  about  1-2x10  cm/sec. 
Typical  ion  drift  velocities  in  atmospheric  air  at  the  outer  limits  of 
the  anode  extension  (E/P  =  25,000/760  -  33)  also  run  in  the  neighborhood 
of  5-8x10^  cm/sec. 

Defining  the  slip  parameter  S  =  u/v.  *  ,  significant  flow 

i 

effects  can  be  expected  when  the  slip  parameter  is  of  order  one,  especial¬ 
ly  in  a  dense  collision  dominated  discharge  where  the  mean  neutral  gas 
motion  is  easily  transfered  to  the  ion  population.  As  can  be  seen  in  the 
speci fi c  example , 


S 


4  A 

u  _  1x10  2x10  _ 

- —T  r  - 

k;E  8x10*  5x10* 


07) 


and  substantial  space  charge  redistribution  should  be  anticipated. 

Going  back  to  Equation  16,  and  rewriting  to  include  flow  effects, 

we  have 

a.  No  electron/ion  coupling: 
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,  (vi+u)  2,e  >  SE('+S)ve  ,e  4 

6u  ’  — 2“^e  (2F5 - 1 - (2F} 


v.v  2  (k.E)v  2 

i  e  ,e  »  _  i  e  ,e 


-(lr) 


<*r> 


O  ;2  2p'  ]2  '2p 

J  J 


08) 


(18a) 


and 


5* 


u  .  ( 1+S) 


'»  3  -i— -  3  ■  i.  .  .m  i 

1  17 


(18b) 


b.  Heavy  electron/ion  coupling  (high  space  charge  fields): 


J  J  ^ 

Ju  Ju 


(19) 


v.3 

5o  •  3-%» 
J 


(1 9a) 


and 


(1+S)' 


(19b) 


whe  re  j  *  *  — 

Jo 


From  these  considerations  it  can  be  seen 
velocity  term  (1+S)  will  vary  with  space  charge 
and  we  will  write  a  generalization  of  1 8b  and  19b 


that  the  power  of  the 
densi ty  (or  current) 
as  follows: 


6* 


( l+S) 


m 


(20) 
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where  m  is  a  variable  exponent. 


Since  continuity  considerations  require  that  the  total  current 
passing  any  equipotential  surface  within  the  discharge  be  constant  the 
total  current  can  be  represented  by  the  electron  current  entering  the 
anode  where  ion  contributions  are  zero,  so 


j 


af {6)6 


(21) 


where  j  is  the  electron  current  flowing  into  the  ionization  region 
(defined  by  x£  which  is  proportional  to  6  )  ,  the  exponential  is 

the  multiplication  factor,  and  f  is  an  as  yet  unspecified  function  of 
<5  ,  therefore  - 


a  fS 

with  flow: 

.  _  .  u  u 

Ju  Joue 

(22a) 

a  fS 

no  flow: 

.  .  o  o 

j  =  j  e 

Jo  Joo 

(22b) 

fa  5  (a*<5*-l) 

then 

•  J.  •  JL  0  0 

J*  **  JQ*e 

(23) 

NOTE:  The  *  value  of  a  quantity  is  used  consistently  to  indicate  the 

value  of  that  quantity  in  a  flow  field  normalized  by  its  no 

flow  value  (except  for  pressure  -  P*  to  be  introduced  later  is 

equal  to  P/P  *  local  value/value  at  STP  =  760  TORR) . 
o 
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The  function  f (5 )  has  been  introduced  to  relate  the  critical 


ionization  distance  (that  point  at  E(x)  -  EAPPUE0  +  ESPACE  CHARGE  Z 

E.)  to  the  anode  extension  6  .  It  is  apparent  that  6  is  proportional 

to  x  i.e. ,  as  <S  increases  either  due  to  convection  or  increasing 
c 

applied  voltage,  x£  also  increases  due  to  increased  applied  or  dis¬ 
placed  space  charge  field,  and  conversely,  as  the  anode  extension  shrinks, 
x^  will  decrease.  However,  due  to  the  non-linearity  of  the  fields 
involved,  a  variation  in  5  will  not  produce  a  linear  change  in  xc 
To  account  for  the  non-linearity,  f (S)  is  assumed  to  obey  a  simple 
power  law  form 

f  =  CS*n  {2k) 

Squaring  j*  and  substituting  into  Equation  20,  yields  the  following: 


6* 


0+sr 


J'o  2exp[(2C6'na06o(a*<S*r.l)] 


(25) 


Rearranging  and  taking  the  natural  log  of  both  sides, 


Hn<5*  +  2Cho<5o5*n(a*<5*-l)  -  m2n(l+S)  -  Hnjo*2 


and  solving  for  o*"1  and  collecting  6*  terms  on  the  left  side, 


s*"  -  ^ *  teV’ik  -  Wr»> <’*s>  -  *"C  <26> 

o  o  o  o 
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This  is  the  basic  anode  extension  equation  in  a  flow  field,  where  again 
the  *  values  are  the  flow  values  compared  to  the  no  flow  values,  i.e.: 

5*  *  5  /S  =  normalized  anode  extension 
u  o 

a*  *  a  /a  =  normalized  ionization  coefficient 
u  o 

S  ■  u/v.  =  the  ion  slip  parameter 

j  *  =  j  /j  =  normalized  saturation  current 

o  ou  oo 

C  =  CONSTANT 

The  exponent  n  of  6*  turns  out  to  be  a  free  parameter  of 
the  problem  whose  exact  nature  is  not  yet  well  established,  but  is  under 
investigation.  The  parameter  m  is  the  power  of  the  velocity  term 
(1+S)  which  analytically  was  determined  to  range  from  approximately  1 
to  3  ,  although  experimental  data  indicate  that  occasionally  a  higher 

value  is  required.  The  reasons  for  this  disparity  are  also  under  active 
investigation. 

Inspite  of  the  limitations  cited  above,  this  equation,  even  in 
its  preliminary  form  (Equation  26)  describes  all  the  previously  observed 
breakdown  voltage  characteristics  on  the  particular  apparatus  used  in 
the  initial  tests  for  various  velocities,  turbulence  levels,  and  flow  ori¬ 
entations  tested,  even  including  some  seemingly  anomalous  points  in 
earlier  tests  which  radically  departed  from  the  norm  as  reported  in 
References  3  and  7.  The  equation  also  predicts  some  salient  features  in 
the  low  velocity  ranges  which  had  not  previously  been  observed  due  to 
strong  emphasis  on  high  speed  testing.  These  low  velocity  features  have 
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been  confirmed  in  recent  tests  and  will  be  described  in  detail  in  part  & 
of  this  section  of  the  report. 


To  check  the  mathematical  correctness  of  the  form  of  this  equa¬ 
tion,  when  there  is  no  flow,  S  «  0  ,  and  by  definition,  the  (*)  terms  5 
1  .  Therefore,  the  three  i,n  terms  go  to  zero,  a*  =  l  and  we  are 
left  with  the  identity  0=0  as  it  should. 

The  peculiar  geometric  factors  of  a  problem  are  all  contained  in 

the  <S  terms.  Gas  and  electric  field  parameters  enter  through  the  a 

and  S  terms;  and  turbulence  effects  enter  through  a*  as  will  be 

explained  in  the  following  sections.  Secondary  electron  production 

processes  are  included  in  j  *  . 

o 

To  obtain  a  solution  to  Equation  26,  requires  some  auxiiliary 
equations  which  are  developed  below. 
k.  Auxiiliary  Equations 

a.  The  no-flow  anode  extension  (5  ): 

o 

The  basic  assunption  of  the  anode  extension  model  is  that 
the  gap  breakdown  potential  is  a  linear  function  of  the  effective  gap 
length,  i.e.,  the  actual  gap  length  less  the  anode  extension  (which  in 
effect  is  a  portion  of  the  gap  already  in  a  state  of  local  breakdown). 
Therefore: 


VbQ  s  d'  z  d  -  So  -  E(d-S0)  +  Vc  (27) 

Vb  =  no  flow  breakdown  voltage 
o 

d"  5  effective  gap 
d  =  actual  gap  length 
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5  =  no  flow  anode  extension 

o 

E  -  mean  electric  field  in  effective  gap 
V£  =  corona  initiation  voltage 

For  a  specific  geometry,  VbQ  can  be  calculated  from  the  condition  that 
E(x)  =  E.  at  x  ■  C-<$0  ,  resulting  in 


VbQ  -  jiP*(R+C-6o)log(4d/R)  (29) 

for  the  hyperbololdal  point-to-plane. 

An  analytical  form  can  also  be  determined  for  Vc  as  in 
Section  2  by  assuming  that  Vc  is  a  linear  function  of  the  point  field, 
i.e.,  when  the  region  in  the  vicinity  of  the  point  reaches  the  reduced 
field  strength  for  ionization.  Again  for  the  hyperboloidal  point,  this 
results  in 


PRlog(4d/R)  Pq 

and 

CiEi 

Vc  a  — -P*Rlog (4d/R) 


(30) 


(31) 


and  inserting  Equations  29  and  31  Into  27  yields: 


5 


o 


d 

P*(£.|-4log(WR)  +  1 


(32) 


b.  Ionization  Coefficient  Equations  (a  and  a*): 

(1)  No  Flow  -  use  Equation  6  with  0  subscript: 


A  fie 


BP 

T“ 


(33) 


(2)  With  Flow/Turbulence  Effects  : 

The  primary  coupling  of  a  with  the  flow  is  result  of 
pressure/density  fluctuations  in  a  highly  turbulent  flow.  Therefore,  by 
substituting  P  ±  AP  into  the  above  equation  and  taking  a  simple  average 
[a(P+AP)  +  a(P-AP)]  11  ,  and  eliminating  second  order  terms  results  in  - 


or 


a 

*  r 


ct*  -  — 
ao 


1 

r 


BAP 

E 


BAP 

E 


+  e 


BAP 

E 


(34) 


+  e 


BAP 

E 


(35) 


The  pressure  fluctuations  due  to  turbulence  can  be  related  to  the  mean 
flow  velocity  through  the  local  pressure-density  relation  in  an  incom¬ 
pressible  flow,  i.e.,  the  total  pressure  equals  the  sum  of  the  static 
and  dynamic  pressures. 
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(36) 


P  -  PD  +  ^  pu2 


Assuming  local  test  section  total  pressure  remains  constant  in  the 
presence  of  turbulence ,  and  that  the  fluctuating  component  of 

velocity,  is  proportional  to  U  ,  the  mean  flow  velocity,  then 


P  -  Pq  ■  iP  :  Cpu' 


The  quantity  (  'J u'2/U)  is  called  the  turbulence  intensity  (Reference 
18).  At  constant  temperature,  p  -  P  ;  therefore  for  constant  T  , 


but  also 


AP  *  SpU 


S  -  U/kjE  U  -  Sk.E 


u/k.E  -*■  u  =  U  ■  Sk.E 


k  *  k  — — 
Ki  Ko  P 


AP  =  EP(Sk.E)' 


Substituting  into  Equation  35  yields  finally  - 


2  2 

B£Pk  s  t 
o  o 


which  when  u  =  U  =  S  =  0  ,  reduces  to  a*  ■  1  as  it  should. 
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c.  Breakdown  Voltage  Equation  with  Flow  and  Turbulence 

The  ultimate  objective  of  the  above  derivations  is  to  arrive 
at  an  expression  for  computing  the  fundamental  gas  dynamic  effects  on  the 
breakdown  voltage  of  an  electric  discharge.  Approximating  the  breakdown 
voltage  ratio  as 


Vb* 


Vb  d-S 
_ u  u 

Vb  '  TT" 
o  o 


then  in  final  form, 


from  Equation  27  ,  (41) 


i  -  4*)s* 

Vb*  -  - - -  (42) 

0  -/) 

which  also  reduces  to  Vb*  *  1  when  U  -  0  <5*  •  1  . 

It  should  also  be  noted  that  when  5*  ■*  0,  Vb*)  -*•  -  .......  ,  the 

max  5  ’ 

'■f 

same  relation  which  was  postulated  in  Reference  3  based  on  the  anode 
extension  concept. 

d.  Summary  of  Equations 

For  convenience  of  reference  in  the  following  discussion, 
the  basic  breakdown  potential  and  auxilliary  equations  are  repeated 
be  low  : 
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(1)  Normalized  Anode  Extension  <5*  ■ 


6*°  •  +  (2 "  (2toV)S^'n(,+S) 

oo  o  o 


Unj  *2 
Jo 


(26) 


(2)  No  Flow  Anode  Extension  dQ  : 


- -g- - 

p*(rr")  ’og(4d/R)  +  i 


(32) 


(3)  No  Flow  Ionization  Coefficient  a  : 


BP 

a  - 

o 


(33) 


(4)  Normalized  Ionization  Coefficient  a*«  — 

ao 


a*  »  y 


2  2  2  2 
B£P  k  VE  -B£P  k  b^E 
o  o  o  o 


T*~ 


p* 


(40) 


Vb 


(5)  Normalized  Breakdown  Voltage  Vb*  - 


Vb*  - 


’  -r5* 

- 3 - 

(I?) 


(41) 
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5.  Qualitative  Discussion  of  the  Equations 

The  key  to  a  clear  understanding  of  the  equations  lies  in  a 
simple  analysis  of  the  normalized  anode  extension  d*  .  The  basic  form 
of  this  equation  is 


&*  - 


m 


In  fact  it  was  derived  from  Equation  20. 


6* 


(1+S)f 

•4.2 


2  2 

But  S  ~  u  and  j*  ~  n*  which  gives 


m 

<5*  ~  L  -  ^ 

*  7? 


m 


which  is  essentially  the  relationship  derived  in  Reference  3.  Thus,  the 
anode  extension  is  a  product  of  two  factors,  the  pure  convective  velocity 
which,  in  normal  flow,  tends  to  lengthen  6*  ,  and  the  inverse  of  the 
charged  particle  density,  which  tends  to  shorten  6*  due  to  faster 
recombination  at  higher  n*  and  consequently  results  in  less  particle 
dri ft. 

In  the  reverse  flow  mode,  the  5*  equation  can  be  approximated 

by: 


6*  -  (1+S)m 


(45) 
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where  u  and  therefore  S  are  negative,  which  decreases  <5*  in  a  rapid 
fashion  with  increasing  negative  velocity.  It  is  interesting  to  note  that 
in  the  reverse  flow  approximation,  a*  cancels  out  and  Vb*  is  not 
influenced  by  turbulence. 

The  breakdown  voltage  behavior  is  opposite  in  sense  to  5*  ,  so 
when  6*  increases,  Vb*  decreases  and  vice-versa.  Note:  frequent 
reference  to  Figure  6  wi 1 1  be  helpful  in  interpreting  the  following.  As 
mentioned  in  the  reverse  flow  mode,  6*  is  quickly  reduced  and  therefore, 
Vb*  immediately  rises  (Part  2  of  Figure  6)  and  quickly  approaches  an 
asymptotic  limit  (Part  1  of  Figure  6)  as  5*  -*•  0  corresponding  to  the 
requirement  to  breakdown  the  entire  physical  gap. 

In  the  normal  flow  mode  (flow  from  anode  to  cathode),  the  initial 
low  gas  velocity  tends  to  lengthen  6*  ,  which  causes  a  decrease  in  break¬ 
down  voltage  (Part  of  Figure  6),  but  due  to  the  exponential  dependence 
of  n*  ,  j*  ,  and  a*  on  U  (through  S)  ,  with  effective  turbulence  the 
a*  terms  in  the  denominator  of  Equation  26  soon  prevail,  causing  5*  to 
decrease  and  Vb*  to  increase  at  higher  velocities  and  turbulence  levels 
(Part  5  of  Figure  6).  The  velocity  at  which  this  transition  occurs  is  a 
direct  function  of  the  effectiveness  of  the  turbulence  introduced.  If  the 
turbulence  is  inhomogeneous  or  ineffective,  the  Vb*  curve  may  not  reach 
a  minimum  in  which  case,  further  increases  in  velocity  may  actually  reduce 
power  input  to  the  discharge  (Part  7  of  Figure  6,  dashed  curve  extension). 

It  is  interesting  to  note  also  that  in  the  normal  flow  mode  with 


effective  turbulence,  after  initially  increasing,  6*  can  again  approach 
0  as  its  absolute  limit  (Part  6  of  Figure  6)  due  to  the  high  recombination 
rate,  thereby  confining  the  normal  flow  Vb*  rise  to  the  same  asymtotic 


limit  as  applies  in  the  reverse  flow  case,  and  again  corresponding  to  the 
requirement  to  breakdown  nearly  the  entire  physical  gap. 

It  should  also  be  noted  that  this  asymptotic  behavior,  which  has 
been  observed  in  many  cases,  (References  2  and  3),  poses  a  fundamental 
limitation  on  the  stabilization  of  the  discharge  potential  by  gas  dynamic 
means.  The  final  breakdown  may  possibly  be  the  result  of  either  of  two 
separate  sequences: 

(a)  the  relatively  uniform  external  "drift"  field  may  reach 
such  a  high  value  that  a  random  electron  avalanche  originating  in  the 
vicinity  of  the  cathode  can  take  enough  energy  from  the  field  to  com¬ 
pletely  cross  the  interelectrode  gap,  resulting  in  an  internal  break¬ 
down  of  the  gas,  or, 

(b)  if  the  cathode  material  has  a  low  work  function  it  may 
begin  emitting  numerous  secondary  electrons  prior  to  the  drift  field 
reaching  the  gas  breakdown  potential  in  which  case  the  breakdown  proba¬ 
bility  will  rise  and  a  cathode  limited  breakdown  may  occur. 

As  turbulence  intensities  are  increased  for  a  given  gap  length, 
the  reverse  flow  region  is  virtually  unaffected,  except  to  hasten  the 
rise  of  Vb*  to  the  asymptote,  while  in  the  normal  flow  direction, 
increasing  turbulence  results  in  higher  a*  values  with  an  eventual 
upswing  in  Vb*  due  to  the  resultant  <5*  reduction.  Figure  7  illus¬ 
trates  the  progressive  modification  of  the  Vb*  curves  with  increasing 
turbulence.  The  asymptotic  behavior  observed  at  high  normal  velocities 
and  turbulence  levels  confirms  the  apparent  rise  in  the  Vb*  curve  and 
explains  logically  the  apparently  anomalous  "short  gap,  high  current" 
data  points  reported  in  References  3  and  7. 
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The  analytic  solutions  illustrated  in  Figure  7  also  provided  an 
additional  area  for  test  and  verification  of  the  breakdown  equations.  It 
can  be  seen  from  the  Figure  that  even  in  the  "short  gap,  high  current" 
cases,  there  should  be  at  least  a  small  region  of  decreasing  Vb*  at  very 
low  normal  flow  velocities  due  to  the  non-turbulent  convective  effect. 
Recent  low  speed  tests  have  confirmed  this  characteristic  as  the  data  in 
Figure  8  show. 

Another  qualitative  assessment  of  the  equations  can  be  made  by 
observing  the  behavior  of  the  curves  as  a  function  of  the  gap  length. 

Figure  9  illustrates  this  behavior  for  a  1  and  5  centimeter  gap  which 
encompasses  the  ranges  tested.  The  trends  of  the  curve  shift  in  both 
normal  and  reverse  flow  with  changing  gap  length  were  as  indicated  by  the 
theoretical  curves. 

Section  6  provides  a  more  quantitative  comparison  of  observed 
data  and  the  appropriate  theoretical  curves.  However,  before  moving  into 
that  subject,  a  very  interesting  and  rather  striking  qualitative  comparison 
can  be  made  of  the  form  of  Equation  26  and  Meek's  Breakdown  Criterion 
(as  presented  in  Reference  12)  for  the  uniform  field  in  a  static  (non¬ 
flowing)  discharge.  As  will  be  seen,  despite  the  very  different  conditions 
in  the  two  discharge  situations,  and  the  very  different  approaches  used  in 
deriving  the  breakdown  criteria,  the  resulting  breakdown  equations  are 
remarkably  simi lar. 

Since  a  description  of  the  non-uniform  field  anode  extension  has 
already  been  developed,  a  short  description  of  the  Meek-breakdown  scenario 
in  the  uni  form  field  will  help  put  the  remai nder  of  this  sect  ion  in 


perspective.  The  Meek  Criterion  is  based  on  the  relative  space  charge 
field  strength  at  the  head  of  an  electron  avalanche  compared  to  the  applied 
electric  field.  If  the  avalanche  travels  a  certain  critical  distance  under 
the  influence  of  the  applied  field,  the  space  charge  field  at  the  surface 
of  the  space  charge  region  will  become  the  same  order  of  magnitude  as  the 
applied  field.  When  the  field  ratio  K  =  E^/E  -►0.1-1,  corresponding  to  the 
avalanche  traveling  the  critical  distance  xc  ,  the  space  charge  field  can 
no  longer  be  ignored.  The  Meek  Criterion  from  Reference  12  is  written  as 
fol lows : 

ax  +  in—  *  14.46  +  in^r^-  -0.55-npx  +  Hn  x  (46) 

P  Pd 

The  method  of  solution  and  some  results  are  discussed  in  detail 
in  the  reference.  What  we  will  concentrate  on  are  some  straight  forward 
substitutions  into  the  above  equation  for  comparison  with  the  anode 
extension  model.  If  the  gap  length  is  adjusted  such  that  an  avalanche 
becomes  critical  just  as  it  reaches  the  anode,  then  the  electric  field 
lines  will  be  distorted  in  the  direction  of  the  space  charge,  as  in 
Figure  10a  and  the  concentration  of  charge  will  produce  a  non-uniform 
field  and,  in  fact,  just  as  the  avalanche  touches  the  anode,  it  acts  as  a 
small  area  (like  a  point)  of  anode  extension.  The  axial  length  of  this 
extension  is  2r  ,  but  r  (the  avalanche  radius)  is  proportional  to  the 
distance  it  has  traveled,  so  from  2r  ~  6  ~  x  , 

x  =  c5  (47) 
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and  it  can  already  be  seen  that  the  situation  looks  very  much  like  the 
point-to-plane  geometry.  After  making  this  substitution,  and  then 
rearranging  and  substituting  non-dimensional i zed  flow  variables  for  all 
the  quantitites,  the  result  is 


a  cx*kp 


-  <« +  (48> 


-  compare  with  Equation  26  for  the  non-uniform  flow  field  - 


«*"  *  Wr»>,s* "  <2 »"V2 


The  transformed  Meek  Criterion  looks  very  much  like  equation  26  with  some 
sign  changes,  but  the  basic  form  is  nearly  identical  to  the  non-uniform 
flow  field  equation.  The  sign  changes  can  be  accounted  for  by  the  obser¬ 
vation  that  in  the  non-uniform  field  equation: 


and  as  the  point  radius  (sphere  or  wire  to  plane  equations  are  similar) 

R  -*•  ®  ,  the  log  term  becomes  negative  and  the  equations  merge  to  nearly 
the  same  form. 

After  seeing  how  similar  the  actual  physical  situations  are  at 
breakdown  though,  the  similarity  should  not  come  as  a  surprise.  In  fact, 
what  the  similarity  further  suggests  is  that  perhaps  there  really  is  no 
absolute  uniform  field  breakdown.  In  a  uniform  gap  at  a  potential  just 
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below  breakdown,  statistically,  numerous  small  electron  avalanches  may 
originate  and  die  out  in  a  short  distance  as  depicted  in  Figure  10a.  Many 
of  these  mini -avalanches  must  meet  the  anode  to  form  minute  deformations 
of  the  field  there,  any  one  of  which,  if  large  enough,  can  ultimately 
trigger  total  gap  breakdown  similar  to  the  non-uniform  point-plane. 

This  brings  us  to  an  interesting  obeservation  concerning  a  highly 
non-uniform  field  situation  where  the  accumulation  of  space  charge  around 
a  positive  point  tends  to  distort  its  apparent  shape  and  radius,  the  result 
of  which  is  to  tend  to  reduce  the  non-uniformity  and  go  toward  a  more  uni¬ 
form  field  configuration.  From  these  two  examples,  it  would  appear  that 
the  "natural"  state  at  breakdown  in  a  discharge  is  a  moderate  state  of 
non-ini formi ty  rather  than  at  either  extreme. 

6.  Quantitative  Comparison  of  Data  and  Theory: 

Figures  11  through  17  present  some  typical  data  recorded  on 
various  runs  with  the  point  and  wi re-to-pl ane  electrodes.  Conditions 
ranged  from  high  reverse  flow  velocities  (-120  m/s)  to  high  normal  flow 
(80-100  m/s  with  turbulence).  No  turbulence  could  be  introduced  in  the 
reverse  flow  mode,  since  as  illustrated  in  Fig.  1,  the  cathode  is  constructed 
of  aluninum  honeycomb  flow  straightener  which  assures  smooth  laminar  reverse 
flow  throughout  the  velocity  range  tested. 

The  data  presented  were  collected  over  a  six  week  period  in 
varying  weather  conditions  and  are  not  corrected  for  relative  humidity  and 
ambient  temperature  effects,  both  of  which  affect  the  properties  being 
measured.  Correction  for  these  effects  might  reduce  the  scatter  evident 
in  various  data  sets,  especially  in  the  composite  summary  plots.  An 
additional  source  of  statistical  scatter  due  to  unequal  current  flow  to 
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various  pins  or  wires  may  also  be  significant,  especially  in  a  discharge 
such  as  ours  which  is  not  ballasted  to  produce  uniform  current  flow  from 
each  element.  As  noted  in  Ref.  3,  however,  effective  turbulence  does 
cause  more  uniform  activation  of  all  pins  and  therefore  reduces  the 
adverse  effects  due  to  no  ballast. 

Figure  11  displays  corona  initiation  and  breakdown  voltages  with 
no  flow  versus  gap  length  for  both  geometries  tested.  Theoretical  curves 
are  also  displayed  which  were  computed  from  the  appropriate  equations 
developed  in  Section  A  .  These  curves  display  an  intersection  of  the 
Vb  and  Vc  characteristics  at  the  shorter  gap  lengths  and  are  very 
similar  in  appearance  to  those  presented  in  Reference  8  for  single 
point-plane  geometries.  The  intersections  define  what  is  termed  the 
"existence"  region  in  Reference  8  for  the  corona  discharge.  To  the 
left  of  this  point,  a  corona  discharge  will  not  form,  and  the  gap  will 
breakdown  (arc)  upon  reaching  the  value  of  Vc  .  To  the  right  of  that 
point  a  corona  will  form  which  enhances  the  discharge  current.  A  common 
empirical  expression  which  has  been  found  to  quite  closely  fit  current 
data  for  various  non-uniform  geometries  can  be  written  as  follows: 

I  -  KCV(V-Vc) 

The  initial  constants  contain  geometric  and  gas  parameters.  This  equation 
highlights  the  dependence  of  I  on  V  and  on  the  difference  between  the 
applied  and  corona  initiation  voltages.  The  existence  regions  for  the 
present  configurations  are  seen  to  be  for  - 

Points-Plane  d  >  0.2  cm 
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d  >  2.0  cm 


Figure  12  presents  a  comparison  of  Vb*)  from  Equation  42 

max 

when  <S*  ■*  0  .  Due  to  power  supply  limitations,  insufficient  points  were 
accessible  in  the  wire/plane  configuration  because  of  higher  breakdown 
voltage  requirements  to  adequately  define  the  wire/plane  curve  so  only  the 
point/plane  data  are  presented.  This  curve  is  interesting  in  that  it 
shows  at  longer  gaps  flow  cannot  affect  breakdown  voltage  as  radically  as 
it  does  in  a  short  gap.  The  general  trends  of  the  Vb*  vs.  velocity 
data  taken  under  similar  conditions  of  humidity  and  temperature  tend  to 
confirm  this.  However  some  contradictions  are  evident  and  this  cannot 
be  considered  a  conclusive  result. 

Figures  13,  14,  and  15  are  comparisons  of  data  and  theory  for 
three  gap  lengths  of  Vb*  vs.  velocity  .  The  values  of  the  parameters 
C,  m,  n,  and  turbulence  index  used  in  the  calculated  curves  are  given.  It 
has  been  observed  that  C  *  1.3  ,  m  *  3  ,  and  n=  8  generally  result  in 
a  good  first  estimate  for  the  basic  curve,  and  then  the  turbulence  index 
(input  of  turbulence  intensity)  is  varied  to  match  the  data.  In  using 
this  technique,  it  has  been  observed  that  long  gaps  usually  require  a 
lower  turbulence  index  even  though  the  actual  turbulence  being  generated 
at  a  given  velocity  is  the  same  as  in  shorter  gaps.  One  possible  explana¬ 
tion  for  this  is  that,  as  pointed  out  by  Reference  2  after  the  turbulencee 
being  generated  is  fully  developed,  it  begins  to  rapidly  decay  with 
time  and  distance,  so  in  a  long  gap,  although  turbulence  levels  entering 
the  discharge  region  are  the  same  as  in  short  gaps,  the  turbulence 
level  leaving  the  discharge  will  have  decayed  more,  resulting  in  a  lower 
average  turbulence  intensity  in  the  longer  gap.  This  would  indicate  that 
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a  weighting  of  the  turbulence  index  with  gap  length  may  be  required. 

Probe  measurements  with  screen  #3  installed  in  the  test  tunnel  show  that 
at  the  highest  attainable  velocities,  the  turbulence  intensity  peaks  just 
at  the  beginning  of  the  discharge  region  (at  the  pin  tips),  which  would 
support  the  above  argument,  and  which  also  may  be  the  key  to  the  success 
of  the  screen  #9  design. 

Figures  16  a,  b  and  17  a  and  b  are  summary  plots  of  Vb  and  lb 
data  vs.  velocity  for  the  points  aid  wires.  While  many  gap  lengths  and 
turbulence  levels  are  included  under  varying  ambient  conditions,  the  plots 
serve  to  point  out  the  main  features  and  general  trends  of  these  parameters 
with  flow  and  turbulence.  The  lines  on  these  curves  are  not  calculated, 
but  are  drawn  in  to  highlight  the  basic  features  of  the  various  data 
collected  and  to  show  the  similarity  with  the  theoretical  curves  dis¬ 
played  elsewhere. 

The  next  section  qualitatively  examines  the  wire/plane  in  more 
detai  1 . 
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S.  Wi re-to-Plane  Geometry 

1 .  Qualitative  Discussion  of  Differences  from  Point-Plane 

One  of  the  significant  changes  in  going  from  the  point-plane 
geometry  to  the  wire-plane  is  a  major  modification  in  the  shape  of  the 
electric  field  strength  distribution  across  a  gap  of  given  dimension  and 
a  fixed  applied  voltage.  While  the  assumed  distribution  (proportional  to 
x  *)  is  still  "hyperbolic"  as  illustrated  in  Figure  3»  at  the  electrode 
surfaces  the  point  field  is  much  higher  and  drops  off  much  more  quickly 
than  a  wire  field  of  the  same  radius  at  the  same  potential.  Since  in 
both  cases  the  cathode  is  grounded,  the  wire  electric  field  in  the  drift 
region  must  be  higher  than  that  of  the  point-plane.  The  wire  radii  used 
in  the  current  research  started  at  0.08  cm,  but  it  was  found  that  to 
operate  within  power  supply  capabilities,  the  wire  size  had  to  be  reduced. 

As  a  result  the  wires  were  incrementally  filed  down  to  a  wedge  shape  of 
approximately  30°  angle  and  maximum  thickness  of  about  0.05  an.  This 
permitted  some  current  enhancement  within  available  power  supply  limits. 

Several  things  are  evident  in  comparing  point  and  wire  distribu¬ 
tions  (see  Fig.  3  )•  First,  for  a  given  gas,  the  ionization  field  strength 

CJr 

E.  determines  the  length  of  the  anode  extension  6  ,  so  i f  E.  is  lower 

i  i2 

than  the  point  of  crossing  (point  P)  of  the  curves  at  a  given  voltage, 

6  will  be  greater  than  S  ,  but  If  E.  is  higher  than  this  point, 
w  p  I  j 

then  6  will  be  greater  than  6  .  Since  as  will  be  discussed  later, 

p  w 

the  current  is  a  function  of  S  ,  this  behavior  can  have  important  con¬ 
sequences. 

Another  important  characteristic  difference  of  the  point  and  wire 
distributions  lies  in  the  transition  region  from  the  high  field  to  the  drift 
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field.  The  transition  is  well  defined  for  a  sharp  point,  but  becomes 
very  slow  and  indefinite  in  wires  and  points  of  large  radii.  This 
results  in  a  less  well  defined  anode  extension  in  the  flatter  distributions, 
which  should  affect  the  flow  interactions  in  the  loosely  bound  outer  por¬ 
tions  of  the  space  charge  (<S)  region.  Changes  in  the  behavior  of  the  <S 
region  with  flow  should  show  up  as  changes  in  Vb*  with  velocity.  Some 
major  changes  in  the  wire-plane  and  point-plane  Vb*  vs.  velocity  curves 
have  been  well  documented  in  the  current  research.  Figure  18  illustrates 
the  basic  nature  of  these  changes. 

It  has  been  consistently  noted  that  the  wire  curves,  while  still 
exhibiting  the  low  velocity  depression  of  Vb*  ,  almost  always  take  an 
early  (i.e. ,  low  velocity)  upturn  toward  the  asymptote  whether  the  flow  is 
turbulent  or  not  (Fig.  8).  The  convective  depression,  caused  by  an  initial 
lengthening  of  6  ,  is  very  small,  and  it  is  thought  that  the  rapid  rise  in 
Vb*  is  due  to  a  rather  large  movement  and  dissipation  of  loose  space 
charge  at  the  outer  limits  of  the  anode  extension.  In  this  region,  due  to 
the  flatness  of  the  E-field,  the  ion  slip  parameter  (S  ■  u/k.E)  can  attain 
a  significant  value  over  much  of  6  at  once,  with  the  result  that  the  outer 
portion  of  the  anode  extension  actually  blows  off  or  becomes  too  tenuous  at 
low  velocity  to  act  as  an  extension,  which  then  causes  a  large  and  rapid  rise 
of  Vb*  .  The  asymptotic  limits  of  Vb*  in  both  the  normal  and  reverse 
flow  directions  are  interpreted  again  as  the  requirement  to  break  down  the 
entire  physical  gap,  as  in  the  point-to-plane  development. 

A  third  characteristic  of  the  wire  geometry  is  the  inherently 
higher  values  of  electric  field  in  the  drift  portion  of  the  gap.  This 
increase  in  the  drift  field  increases  the  probability  that  high  ionization 
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and  streamer  activity  in  the  region  near  the  wire  will  be  able  to  develop 
and  sustain  self-propagation  through  the  drift  field  region  to  the  cathode 
by  the  streamer  mechanism  discussed  in  most  texts  on  the  subject  (see  Refer¬ 
ences  8,  12,  15,  16)  resulting  in  breakdown.  On  the  other  hand,  however,  in 
a  discharge  designed  for  optimum  energy  transfer  to  a  laser  gas  (pumping), 
the  design  of  the  electrodes  will  to  some  extent  be  dictated  by  this  pumping 
requirement.  In  practice,  a  trade-off  is  required  between  electrode  design 
optimized  for  total  power  input  to  the  gas  and  that  design  optimized  for 
best  E/N  in  the  drift  region  for  most  efficient  pumping.  In  an  electric 
discharge  laser  cavity,  for  instance,  an  input  power  of  1000  watts/cm^ 
means  very  little  if  only  two  percent  is  going  into  the  lasing  states  due 
to  an  inefficient  E/N  in  the  pumping  region. 

One  of  the  attractive  features  of  the  wi re-to-pl ane  geometry  is 
the  increased  ionization  volume  available  continuously  along  the  wire.  The 
point- to- pi ane  ionization  volume,  on  the  other  hand,  is  confined  to  what 
would  amount  to  discrete  points  along  the  wire. 

The  point-plane  produces  discrete  3"dimensiona1  axially  symmetric 
multiplication  volumes,  at  each  point,  while  the  wire-plane  system  is  a  2- 
dimensional  continuous  volume  spanning  the  length  of  the  wires.  Most  elec¬ 
tron  multiplication  occurs  within  these  volumes,  and  so  it  seems  as  i f  a 
wire  should  produce  more  current.  But  the  problem  is  complicated  by  the 
electric  field  distribution  which  determines  the  longitudinal  dimension  and 
a  the  multiplication  coefficient.  As  will  be  seen  total  current  is  not 
a  function  of  a  and  <5  alone,  but  a  function  of  the  product  (os6)  or 
(a *6*)  in  the  flow  situation. 

2.  Total  Current  Considerations 

As  pointed  out  in  Part  II,  Section  3,  since  the  total  current  at 
the  anode  surface  is  due  to  electrons  entering  the  anode  from  the 
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multiplication  region  (the  positive  ions  being  repelled),  the  total  discharge 
current  can  be  expressed  as  a  function  of  the  number  of  electrons  entering 
and  their  drift  velocity  (the  assumption  is  made  that  if  negative  ions  are 
present,  they  give  up  their  electron,  soon  after  drifting  into  the  high 
field  region;  this  assumption  is  supported  by  Reference  15  and  some 
sample  calculations  presented  in  Reference  8). 

A  single  electron  entering  the  high  field  region  will  exponentially 
multiply  to  produce 

f°  adx 
x 

n  ■  n  e 
o 

electrons  at  the  anode  surface  with  nQ  *  1  .  Since  the  current  density  is 
the  current  per  unit  area, 

Je  “  cneve  *  i/A  (50) 

the  total  current  (which  must  be  constant  across  the  discharge)  can  be 
expressed  as 

_  f°  adx 

'  ■  1  Qe  “  jeA  (50 

This  allows  the  total  current  to  be  discussed  in  terms  of  electrons 
arriving  at  the  anode.  That  number  is  strongly  dependent  upon  how  many 
electrons  enter  the  multiplication  region,  whether  as  free  electrons,  or 
attached  to  atoms  or  molecules  as  negative  ions.  The  electron  input  is 
described  as  iQ  in  the  above  equation  and  is 

i  ■  (e  n  v  +  e  n  v  )A  =  ye  (52) 

o  e_  e 
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where 


n  =  free  electron  density  from  drift  region 
o 

n_  =  negative  ion  density  from  drift  region 
v  ,v.  =  particle  drift  velocities 

C  I 

A  =  area  through  which  particles  pass  into  multiplication  region 
Y  =  secondary  electron  supply  rate. 


Without  covering  many  details  or  methods  already  discussed  in 
Section  A,  the  aoove  considerations  lead  to  the  normalized  current  equation 
below 


i  * 


i 
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Y*e 


2Ca  5  S*n(a*$*-1) 
o  o 


(53) 


where 

Y*  =  i  /'  =  normalized  secondary  electron  supply  rate. 

From  equation  52  it  can  be  seen  that  the  units  of  y  are  electrons/second, 
and  that  it  is  an  electron  flow  rate  into  the  multiplication  region.  It  is 
the  counterpart  of  the  Townsend  secondary  ionization  coefficient  in  a  non¬ 
flowing  discharge. 

Understanding  the  nature  of  y  »  are  ^  able  to  analyze 
possible  flow  effects.  Electrons  coming  into  the  ionization  region  can 
arrive  in  two  forms,  either  as  free  electrons  at  the  characteristic 
electron  drift  velocity  which  is  generally  much  higher  than  typical 
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flow  velocities,  or  attached  to  a  heavier  particle  in  the  form  of  negative 
ions,  etc.,  whose  drift  velocities  and  therefore  delivery  rates  can  be 
greatly  altered  by  the  flow. 

Now  we  discuss  the  third  component  of  current,  the  positive  ion 
contribution.  In  the  determination  of  total  current  by  counting  electrons 
flying  past  the  anode  surface,  the  positive  ions  play  an  indirect  role  by 
supplying  secondary  electrons  primarily  through  impact  with  the  cathode  where 
energetic  collions  can  result  not  only  in  neutralization  of  the  ions,  but 
also  release  of  additional  electrons  which  then  begin  to  move  toward  the 
anode.  These  secondary  electrons  may  proceed  as  free  electrons  or  attach 
to  particles  and  proceed  as  negative  ions,  thus  the  positive  ion  contri¬ 
bution  shows  up  as  an  enhancement  of  both  and  yn„  •  Since  both 

positive  and  negative  ion  drift  velocities  are  of  the  same  order  of  magnitude 
as  the  flows  being  discussed,  it  can  be  anticipated  that  the  total  current 
would  be  significantly  affected  by  flow  through  ion  slip  -  and  electron 
trapping  discussed  in  Part  I  I  A. 

Considering  the  free  ele^-t^on  component  alone,  the  total  input  of 
electrons  to  the  multiplication  zone  is,  from  equation  52 

i  ■  (e  n  v  )A  (54) 

o  e  e 

At  first  it  would  seem  that  electron  velocity  in  the  drift  region  should 
not  be  significantly  affected  by  the  flow  due  to  low  electron  slip  and 
negligible  trapping  due  to  space  charges.  However 

vft  s  k^E (x)  (55) 
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In  nearly  all  non-uniform  geometries,  the  electric  field  is  given  by 
expression  of  the  following  form  - 


E(x) 


-  i  cV 
7  log ( CONST) 


with  the  result  that 


E*(x)  *  V* 


and  at  breakdown, 


Eb*  =  Vb* 

Therefore  from  equations  $k,  55,  and  58  - 
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This  equation  shows  that  even  free  electron  velocities  can  be  altered 
indirectly  by  flow  modifications  which  affect  the  breakdown  voltage 
and  hence  the  electric  field  even  in  regions  where  electron  trapping 
not  a  significant  factor.  Vb*  and  therefore  electron  speed  can  be 
by  }Q%  or  more  at  moderate  flow  velocities. 

Now  the  total  current  can  be  written  as  follows  - 
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The  n  *  term  is  a  function  of  the  total  current,  and  all  associated 
e 

secondary  electron  supply  processes  such  as  cathode  emission,  photoioniza¬ 
tion,  etc.,  whose  variations  with  flow  are  as  yet  undetermined. 

In  a  discharge  where  many  of  the  secondary  electrons  are  liberated 

by  positive  ion  impact  with  the  cathode,  however,  n  *  will  be  a  function 

e 

of  the  ion  slip  and  the  current  equation  may  assume  the  following  form  - 


i*  -  Vb* 


(1+S)  +  y* 
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(61) 


Neglecting  y*  temporarily,  Equation  61  will  be  used  to  qualita¬ 
tively  examine  the  behavior  of  the  breakdown  current  with  flow  and  turbu¬ 
lence  and  to  attempt  to  explain  some  of  the  observed  differences  in  the 
point-to-plane  and  wi re-to-plane  Vb*  and  lb  vs.  velocity  data  collec¬ 
ted  during  the  current  research  effort. 

Figures  16,  17  and  18  are  typical  plots  of  breakdown  voltage  and 
current  vs.  velocity  for  the  point  and  wire  configurations  and  will  be 
referenced  frequently  in  the  following  discussion, 
a.  Point-Plane  Current  Interpretations 

Since  breakdown  voltage  variations  with  flow  were  already 
discussed  in  detail  in  Part  II. A. 5,  we  will  concentrate  on  some  typical 
point-plane  current  data  first.  Beginning  at  the  origin  (U-0)  in  Fig. 
17a,  it  can  be  seen  that,  as  velocity  increases  in  the  normal  flow 
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(positive)  direction,  the  current  begins  to  rise  (region  1),  reaches  a  low 
maximum  (region  2)  and  then  declines  slightly  (region  3).  A  typical 
Vb*  curve  (Fig.  8)  shows  a  depressed  breakdown  voltage  in  this  region, 
which  from  Equation  61  tends  to  decrease  lb  but  the  increasing  velocity 
term  helps  to  compensate  for  the  reduction  as  does  the  product  (a*<$*)  in 
the  exponential,  and  the  current  in  regions  1  thru  3  is  dominated  by  the 
behavior  of  5*  (in  this  velocity  region  turbulence  intensity  is  not  yet 
sufficient  to  cause  significant  enhancement  of  a*  ).  As  velocity  increases, 
if  an  effective  turbulence  generator  is  in  place,  a*  becomes  dominant  due 
to  its  rapid  exponential  growth  with  increasing  velocity  (Equation  40) ,  and 
in  accordance  with  Equation  61  an  overexponential  rise  in  lb  is  noted 
(region  4  of  Fig.  17a).  The  figure  also  shows  Vb*  leveling  at  the  normal 
flow  asymptote  which  would  also  cause  a  leveling  of  the  lb  character  st i c 
or  at  least  a  leaning  toward  a  (1+S)  linear  rise  in  accordance  with  Equa¬ 
tion  61. 

b.  Reverse  Flow  Current 

Interpretation  of  the  reverse  flow  (negative  velocity)  current 
and  voltage  characteristics  in  the  point-plane  and  wire-plane  is  essentially 
the  same,  so  the  following  discussion  applies  to  both  configurations.  Since 
no  turbulence  was  introduced  in  the  reverse  flow  mode,  other  processes  tend 
to  dominate  breakdown  characteristics  in  this  regime.  Equation  1*5  shows  a 
rapid  decrease  in  6*  with  reverse  flow  which  causes  an  equally  rapid  rise 
in  Vb*  ;  the  Vb*  increase  and  (1  +  S)  velocity  term  decrease  might  be 
expected  to  neutralize  each  other  as  in  the  normal  flow  mode,  leaving  the 
major  influence  on  lb  the  6*  term  again.  Since  5*  is  rapidly  decreas¬ 
ing,  this  line  of  reasoning  would  have  lb  decreasing  at  a  fairly  rapid 
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race.  However,  Che  experimenCal  daCa  consisCendy  show  a  sCeep  rise 
in  lb  aC  relacively  low  negacive  velociCies  (region  5  of  Fig.  17a),  and 
a  leveling  or  decrease  aC  higher  velociCies  (regions  6  and  7).  The  high 
velocity  reducCion  in  lb  is  easily  underscood  via  Equacion  61,  since  the 
velocity  term  (1+S)  will  eventually  approach  zero  at  high  enough  nega¬ 
tive  velocities  causing  lb  to  approach  zero.  Also  since  S*  is  con¬ 
stantly  decreasing  with  increasing  reverse  flow,  this  tends  to  reduce  lb 
more  at  high  velocity.  (In  one  experimental  case  lb  was  reduced  to 
approximately  1/10  of  its  no  flow  value).  It  seems  that  the  only  factor 
going  in  favor  of  increasing  reverse  flow  current  in  the  low  velocity 
region  is  the  increase  in  Vb*  ;  but  numerical  calculations  show  that  the 
velocity  term  more  than  offsets  this  factor.  Careful  review  of  the  equa¬ 
tions  brings  aid  from  an  unexpected  source  to  explain  this  strong  low 
negative  velocity  current  enhancement  through  the  a*  term. 

Recalling  the  form  of  a  from  Equation  6,  it  can  be  seen  that 
a  is  a  function  of  E  as  well  as  P  ,  and  as  Vb*  increases  in  the 
reverse  flow  mode,  E  also  increases,  since  from  Equation  1  E  -  V  . 
Mathmati ca 1 ly  this  leads  to  the  following  enhancement  of  a  due  to  the 
rise  of  E  for  the  point-plane  - 

a*(AE)  2  exp£(-BP/Vbo){(R  log(4d/R)l(l/Vb*-l) }  (62) 

Another  subtle  effect  may  also  come  into  play  here,  since,  as  flow  is 
directed  into  the  quasi-steady  no-flow  anode  extension  region,  the  initial 
hydrodynamic  effect  will  be  to  promote  a  flattening  and  radial  spreading 
of  the  dense  space  charge  boundary.  In  this  case  the  point  radius  R 
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should  be  replaced  by  ^£pp  which  will  be  a  complicated  function  of  the 
flow  and  electric  fields.  Qualitatively,  however,  it  is  expected  that  after 
initial  low  velocity  spreading  which  enhances  a*  (via  Equation  62)  and 
therefore  i*  ,  i*  will  eventually  decline  due  to  the  continued  decreases 
in  6*  and  the  (1+S)  velocity  term  in  accordance  with  Equation  61  and 
R^PP  will  also  decline.  The  end  result  of  these  processes  is  a  rapid 
reverse  flow  current  rise  and  fall  as  illustrated  by  Fig.  7a.  These  two 
reverse  flow  effects  a*(AE)  and  R  -*■  R^pp  are  not  peculiar  to  reverse 
flow,  but  in  the  normal  mode  they  must  compete  with  other  strong  influences 
which  may  tend  to  mask  their  effects  except  under  special  circumstances. 

For  example,  in  the  normal  mode,  gas  flow  and  positive  ion  repulsion  from 
the  anode  act  in  the  same  direction  on  the  high  density  positive  space 
charge,  so  gas  convection  by  itself  does  not  promote  radial  spreading  as 
it  does  in  the  reverse  flow  mode  where  electrostatic  repulsion  and  convec¬ 
tion  oppose  each  other.  If,  however,  the  normal  flow  becomes  highly  tur¬ 
bulent,  the  random  gas  velocity  and  density  fluctuations  in  the  ionization 
region  tend  to  enhance  ionization  and  the  transport  properties  of  the  gas 
resulting  in  faster  radial  diffusion  and  dynamic  spreading  of  the  positive 
charges,  and  as  in  the  reverse  flow  case,  the  effective  radius  of  the 
virtual  anode  is  increased.  Reference  3  presented  photographs  of  various 
discharge  conditions  which  visually  illustrates  this  turbulent  enhancement, 
the  glow  region  in  the  turbulent  case  being  3"5  times  the  radius  of  the 
no-flow  or  laminar  flow  cases. 

c.  Wire-Plane  Normal  Flow  Breakdown  Current  and  Voltage 

Breakdown  data  for  this  configuration  was  limited  to  a  short 
gap  length  range  above  the  existence  point  for  the  corona  discharge  (-2  cm) 
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due  to  power  supply  problems  which  limited  available  voltage  to  approxi¬ 
mately  33 KV.  This  limitation  effectively  restricted  wire-plane  gap  length 
to  2.8  cm  or  less  with  flow.  Since  this  is  less  than  1  cm  into  the 
corona  existence  region,  it  is  not  possible  to  give  a  full  accounting  of 
flow  effects.  In  general,  however,  non-turbulent  curves  were  usually 
flat  with  little  enhancement  of  current,  while  most  curves  with  turbulence 
screen  #9  in  place  showed  slow  but  steady  improvement  with  velocity.  One 
interesting  test  was  performed  with  a  gap  length  below  the  corona  existence 
point,  and  the  introduction  of  turbulent  flow  resulted  in  the  flow  of  a 
small  current  before  breakdown.  This  is  most  likely  due  to  the  pressure/ 
density  fluctuations  in  the  vicinity  of  the  wires  which  cause  corresponding 
fluctuations  in  E/P  .  If  the  lowest  pressure  variations  result  in  a  high 
enough  value  of  the  reduced  field  (E/P),  some  ionization  will  begin  to 
occur  before  breakdown. 

The  wire-plane  breakdown  voltage  curves  were  mentioned  in  Section  II 
A. 6  and  the  general  trends  shown  in  Figs.  8  and  18.  One  additional  obser¬ 
vation  will  be  made  here.  It  has  been  noted  in  the  point-plane  configura¬ 
tion  in  the  normal  flow  mode,  that  as  gap  length  is  shortened  toward  the 
corona  existence  point,  breakdown  voltage  will  rise  at  lower  flow  veloci¬ 
ties.  In  the  most  extreme  point-plane  case,  Vb*  had  reached  the  normal 
flow  asymptote  at  40  m/s  with  no.  turbulence.  In  fact  the  curve  looked 
almost  like  a  wire-plane  curve.  This  seems  to  infer  that  if  gap  length  is 
in  the  corona  existence  region,  but  close  to  the  existence  point,  the 
normal  flow  Vb *  characteristic  may  tend  toward  the  wire-plane  shape  of 
Fig.  18.  That  is,  if  the  gap  is  short  relative  to  the  corona  existence 
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point  a  low  velocity  rise  in  Vb  will  occur.  Additional  data  at  various 
gap  lengths  are  necessary  to  confirm  this  preliminary  observation.  if, 
in  fact,  this  is  a  true  characteristic,  it  is  most  likely  related  to  the 
variation  of  the  axial  breakdown  electric  field  distribution  as  the  gap 
length  is  reduced. 

The  remaining  sections  review  requirements  for  the  optimum 
pumping  of  a  laser  gas  mixture  vs.  electrode  geometry,  a  final  over¬ 
view  of  data  and  theory  developed  earlier,  and  some  conclusions  and 
recommendations  for  further  research. 


C.  Electrode  Design/Optimization 


In  electric  discharges  designed  for  laser  pumping,  two  basic  require¬ 
ments  must  be  satisfied.  First,  the  electrodes  must  efficiently  produce 
a  sufficient  electron  population  to  permit  effective  energy  transfer  to 
the  laser  gas.  Secondly,  the  electrode  design  must  be  such  that  the 
reduced  electric  field  (E/N  or  E/P)  in  the  pumping  region  optimizes  elec¬ 
tron  energy  for  maximum  energy  input  to  the  desired  vibrational  or 
electronic  energy  levels.  Degradation  in  either  of  these  critical  areas 
will  result  in  reduced  laser  power  output. 

Discharge  cavity  size  limitations  or  specifications  will  most  likely 
dictate  the  maximum  gap  spacing  between  anode  and  cathode,  so  within  that 
limitation  the  combination  of  the  stressed  electrode  radius  and  gap  length 
must  be  determined  which  best  satisfies  the  requirements  above.  Since  the 
E-field  distribution  is  a  function  of  both  parameters,  they  cannot  be 
considered  separately.  In  addition,  it  has  been  shown  that  space  charge 
effects  via  the  “anode  extension"  and  radial  spreading  cause  significant 
changes  in  the  “apparent"  electrode  geometry  which  should  be  taken  into 
account  in  the  final  design.  The  gasdynamic  interactions  with  these  space 
charge  induced  dimensional  changes,  in  effect,  provides  an  extra  degree  of 
freedom  in  a  fixed  geometry  cavity  through  some  control  of  the  radial 
spreading  and  anode  extension.  This  in  turn  permits,  to  some  extent,  ad¬ 
justment  of  E/N  and  electron  production  via  flow  and  turbulence. 
Additional  benefits  arise  through  the  enhancement  of  transport  properties, 
improving  stability  while  concurrently  increasing  electron  production. 

Calculation  of  the  reduced  fields  via  the  adaptation  of  the  hyper- 
boloidal  point-plane,  and  an  assumed  wire-plane  distribution  indicates 
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E/N  values  in  the  drift  regions  of  the  discharges  tested  slightly  higher 
than  the  optimum  range  specified  in  Reference  i,  a  condition  which  can  be 
remedied  by  proper  adjustment  of  the  geometric  parameters,  provided  that 
an  adequate  power  supply  is  available. 


III.  SUMMARY 


A.  Data 

The  Vb*  data  for  both  configurations  are  quite  similar  in  the 
reverse  flow  modes,  showing  a  fast  rise  to  an  apparent  asymptotic  limit 
(Figs.  16  a,b).  This  is  presumed  due  to  collapse  of  the  anode  extension 
by  the  gas  flow. 

The  normal  flow  mode  shows  Vb*  initially  decreasing  for  all  gaps 
and  geometries  with  an  eventual  upswing  which  depends  strongly  on  turbu¬ 
lence  intensity,  and  appears  to  also  depend  on  gap  length  relative  to  the 
corona  existence  point,  with  early  (low  velocity)  rise  strongly  favored  in 
the  high  turbulence  and/or  short  gap  cases  (Figs.  8,  16  a,  b,  1 8) . 

The  breakdown  current  data  in  the  reverse  flow  illustrated  similar 
behavior  in  the  point  and  wire-plane  configurations,  starting  with  a  sub¬ 
stantial  low  velocity  rise  followed  by  a  decline  at  higher  velocities. 

The  point-plane  normal  flow  current  behavior  exhibits  several  features 
under  various  conditions.  The  most  promient  of  these  are  outlined  in 
Fig.  17a.  Many  smaller  scale  perturbations  are  evident  in  some  individual 
sets  of  data,  but  the  high  turbulence,  higher  current  curves  seem  to 
follow  the  same  basic  pattern  in  the  point-plane  configuration.  The  wire- 
plane  current  curves  exhibit  much  more  scatter  and  trends  are  not  as  well 
defined  as  with  the  point  geometry.  This  is  probably  attributable  to  the 
fact  that  data  were  not  accessible  in  a  wide  range  of  gap  lengths,  but 
were  restricted  to  gaps  close  to  corona  initiation  where  Vb  and  V c  are 
not  far  apart,  and  therefore  current  levels  are  relatively  low. 
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B.  Theory 


The  theory  developed  in  Part  1 1  accounts  for  the  major  features  of 
the  Vb  and  lb  curves  described  above  as  can  be  seen  in  the  various 
figures,  but  many  of  the  analytical  terms  can  be  improved  and  it  is 
believed  that  further  study  will  lead  to  analytic  expressions  for  the 
various  parameters  of  the  theory  which  have  been  mentioned.  The  clear 
physical  basis  of  the  theory  lends  to  easy  interpretation  of  various  flow 
related  phenomena  and  should  prove  useful  in  visualization  and  design  of 
flow-type  systems  and  possibly  in  trouble-shooting  of  existing  flow  system 
p  rob 1 ems . 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

1.  The  basic  anode  extension  model  developed  herein  appears  to 
offer  a  simple  approach  to  the  interpretation  and  calculation  of  basic 
gasdynamic  effects  on  high  pressure  corona  discharges. 

2.  The  flow  field  interactions  discussed  are  basic,  so  the  approach 
used  in  modelling  the  positive  point-plane  electrode  system  should  be 
easily  applied  to  other  gap  configurations,  polarities,  gas  mixtures,  and 
pressure  regimes. 

3.  A  more  rigorous  evaluation  of  various  terms  in  the  basic  model 
should  provide  better  insight  into  the  dominant  forces  in  various  flow 
regimes,  and  most  likely  will  result  in  definition  of  some  of  the  small 
scale  features  not  yet  well  explained. 

B.  Recommendations 

1.  Continue  data  collection  on  the  wire-plane  electrode  configura¬ 
tion  when  a  larger  power  supply  is  available  (50KV,  50ma  should  be 
sufficient  with  present  hardware). 

2.  Develop  a  good  analytical  expression  for  the  effective  electrode 
radius  for  insertion  into  the  present  point-plane  computational  model. 

3.  Include  a*(AE)  and  radial  spreading  effects  in  normal  flow 
computational  scheme. 

4.  Test  point-plane  in  reverse  flow  mode  with  turbulence  at  low 
velocities  to  verify  Equation  45. 

5.  Develop  comprehensive  analysis  of  flow  on  y*  for  various  types 
of  gas  mixtures,  i.e.,  attaching,  non-attaching,  etc. 
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Figure  1.  Electrodes  and  Turbulence  Generating  Screen  #9. 
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Figure  4.  ELECTRIC  FIELD  AND  VOLTAGE  DISTRI 
BUTION  BETWEEN  PARALLEL  PLATES 
Subscripts :  ( I )  No  Space  Charge 

(2)  Moderate  Space  Charge 

(3)  High  Space  Charge 


A  NOPE  PINS 

ROW  1 
(a) 


ROW  2 
(b) 


ROW  3 
(c) 


Figure  5. 


1 .  Applied  Field 

2.  Space  Charge  Field 

3.  Total  Field 


(o )  Schematic  of  Field  Lines  and  Eguipotentials. 

(b)  Apparent  Dimensional  Changes  Due  to  Space 
Charge . 

(c)  Applied,  Space  Charge  and  Total  Field 
Distributions. 

XQ  and  Xc  are  Critical  Distances  From 

Anode  Where  Electron  Avalanching  Begins 
(E  =  E.  )  With  and  Without  Space  Charge 

Effects. 
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Figure  6.  BREAKDOWN  VOLTAGE  RATIO 
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Figure  7.  PLOT  OF  BREAKDOWN  VOLTAGE  VERSUS  FLOW  VELOCITY 
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Figure  9.  PLOT  OF  BREAKDOWN  VOLTAGE  VERSUS  FLOW  VELOCITY. 


Figure  10.  { o )  Activity  in  a 

(b)  Modified  E 


k. 


FOR  POINT- PLANE  GEOMETRY. 


U  (m/s ) 

PLOT  OF  BREAKDOWN  VOLTAGE  VERSUS  FLOW  VELOCITY 


68 


Vb* 


U(m/s ) 

BREAKDOWN  VOLTAGE  VERSUS  FLOW  VELOCITY 


Figure  16a.  POINTS  TO  PLANE  Vb*  Vs.  U  DATA  SUMMARY. 


TURB  NO  TURB 


Figure  17a.  BREAKDOWN  CURRENT  Vs.  FLOW  VELOCITY  POINT-  PLANE 
SUMMARY  PLOT. 
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